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Abstract 

Macroscopic derivation of the entrainment in in a Andreev-Lifshitz Supersolid induced by 
a surface elastic wave propagating along the flexible interface is conducted by considering the 
nonlinear coupling between the interface and the rarefaction effect. We obtain the critical bounds 
for zero-volume-flow-rate states corresponding to specific rarefaction measure and wave number 
which is relevant to the rather small critical velocity of supersolid flows reported by Kim and Chan. 
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In 1969, it was conjectured by Andreev and Lifshitz 1 that at zero temperature, delocal- 
ized defects may exist in a quantum solid, as a result of which the number of sites of an 
ideal crystal lattice may not coincide with the total number of particles. Originally, this 
conjecture was proposed for three dimensional quantum solids made of atoms ( 3 He, 4 He, 
■ • ■) which do not interact via Coulomb repulsion. The proposed supersolid phase is believed 
to occur due to the quantum behavior of point defects, namely vacancies and interstitials, in 
this crystal of bosons 2-3 . Researchers have found that a small lattice model does not exhibit 
the mesoscopic signature of an intermediate phase separating the solid from the liquid, where 
the solid and the fluid would coexist 4 . Such a vacancy-solid phase was indeed suggested 1 by 
Andreev and Lifshitz if the zero point motions of certain defects become sufficient to form 
waves propagating inside the solid. 

Castaing and Nozieres have later considered 5 such a possibility for spin polarized 3 He. The 
statistics of the defects depend on their nature. For simple vacancies in the crystal, their 
statistics is given by the statistics of the particles out of which the solid is made. If the 
defects are bosons, they may form a condensate, giving rise to a superfluid coexisting with 
the solid. This supersolid phase is discussed in certain bosonic models 6 . If the defects are 
fermions, they may form a Fermi liquid 7 coexisting with the solid, such that the system 
is neither a solid, nor a liquid. Two kinds of motion are possible in it; one possesses the 
properties of motion in an elastic solid, the second possesses the properties of motion in a 
liquid. This interesting issue motivates our present study. 

Early theoretical work by Andreev and Lifshitz 1 and Chester 2 showed that solids may feature 
a Bose-Einstein condensate of vacancies (or interstitial atoms) and thus possess superfluid 
(SF) properties. Quite recently one description of the quantum solid is as a density wave 
that has formed in the quantum fluid 8-9 . The periodicity of this density wave need not 
match precisely to the particle density, so that the ground state may be incommensurate, 
with unequal densities of vacancies and interstitials. Whether or not the same is true for 
quantum fluctuations is not clear at this point. We noticed that previous theories imply a 
corresponding vacancy contribution to the specific heat that is as large as the phonon contri- 
bution near 1 Kelvin 3,10 . Based on these considerations or phenomenological approaches, in 
this letter, we shall demonstrate that wavy flexible interface (between atoms and free vacan- 
cies or defects) or highly-pressured environments 8 can produce elastically deformed interface 
or peristaltic motion will induce time-averaged transport in a Andreev-Lifshitz supersolid. 
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Theoretical studies of interphase nonlocal transport phenomena which appear as a result of 
a different type of nonequilibrium representing propagation of a surface elastic wave have 
been performed before 11-12 . These are relevant to particles flowing along deformable elastic 
slabs with the dominated parameter being the Knudsen number (Kn = mean-free-path/L^, 
mean-free-path (mfp) is the mean free path of the particles, Ld is proportional to the dis- 
tance between two slabs) 13-15 . The role of the Knudsen number is similar to that of the 
Navier slip parameter N s (— fiS/Ld] S is a proportionality constant as u s = St, r : the shear 
stress of the bulk velocity; u s : the dimensional slip velocity; for a no-slip case, S = 0, but 
for a no-stress condition. S — oo, /i is the viscosity). 

We shall choose a periodic domain to simplify our mathematical treatments. The flat inter- 
face is presumed. We adopt the macroscopic or hydrodynamical approach and simplify the 
original system of equations (related to the momentum and mass transport) to one single 
higher-order quasi-linear partial differential equation in terms of the unknown stream func- 
tion. In this study, as the temperature is rather low and the phase is related to the supersolid 
(there might be weakly friction or shearing dissipation in-between) we shall assume that the 
governing equations are the incompressible Navier-Stokes equations which will be associated 
with the microscopically slip velocity boundary conditions along the interfaces 13-15 . To con- 
sider the originally quiescent environment for simplicity, due to the difficulty in solving a 
fourth-order quasi-linear complex ordinary differential equation (when the wavy boundary 
condition are imposed), we can finally get an analytically perturbed solution and calculate 
those physical quantities, like, time-averaged transport or entrainment, critical forcing cor- 
responding to the freezed or zero-volume-flow-rate states. 

We consider a two-dimensional region of uniform thickness. The flat-plane boundaries or 
interfaces are rather flexible and presumed to be elastic, on which are imposed traveling 
sinusoidal waves of small amplitude a (possibly due to quantum fluctuations). The vertical 
displacements of the upper and lower interfaces (y = Ld and — Ld) are thus presumed to be 
i] and — T), respectively, where 77 = acos[27r(x — ct)/X], A is the wave length, and c the wave 
speed, x and y are Cartesian coordinates, with x measured in the direction of wave propa- 
gation and y measured in the direction normal to the mean position of the interfaces. We 
have a characteristic velocity c and three characteristic lengths a, A, and Ld- The following 
variables based on c and L d could thus be introduced : x' = x/L d , y' = y/L d , u' = u/c, 
v' = v/c, rj' = rj/Ld, ip' = ip/(cLd), t' = ct/h, p' = p/(pc 2 ), where ip is the dimensional 
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stream function, u and v are the velocities along the x- and y-directions; p is the density, p 
is the pressure. The primes could be dropped in the following. The amplitude ratio e, the 
wave number a, and the Reynolds number Re (representing the weakly friction or shearing 
dissipation using a viscosity v) are defined by e = a/La, a = 27rL^/A, Re = cL&jv. We shall 
seek a solution in the form of a series in the small parameter e : ip = ipo + eipi + e 2 ip2 + • • 
dp/dx = (dp/dx) + e(dp/dx)i + e 2 (dp/dx) 2 + ••; with u = dip/dy, v = —dip/dx. The 2D 
(x- and y-) momentum equations and the equation of continuity 16-17 could be in terms of 
the stream function ip if the p-term is eliminated. The final governing equation is 

and subscripts indicate the partial differentiation. Thus, we have 

-|v 2 ^ + ^oyVVox - ^OxVVoy = ^ V Vo, ( 2 ) 

d 1 

— VVi + VfyVVis + V'lyV 2 ^ - iioxV 2 ip ly - ^i x V 2 ^ y = ^VVi> (3) 

and other higher order terms. Microscopic boundary conditions imposed by the symmetric 
motion of the interfaces and the non-zero slip velocities 13 " 14 are : u — Kn du/dy, v = 
±dr)/dt at y = ±(1 + rj), here Kn=mfp/L d . The boundary conditions are expanded in 
powers of rj and then e. These equations, together with the condition of symmetry and 
a uniform (dp/dx) , yield : ip = K [(l + 2Kn)y — y 3 /3], K = Re(—dp/dx) /2, *0i = 
{0(y)e ia ^ _ ^ + <j)*(y)e~ ia ( x ~ t " > }/2 , where the asterisk denotes the complex conjugate. A 
substitution of ip\ into Eqn. (3) yields 

d 2 d 2 
{— -a 2 + iaRe\l - K (l - y 2 + 2Kn)l}(— - a 2 )(j) - 2iaK Re = 
dy 2 dy 2 

or if originally the environment is quiescent : K = (this corresponds to a free pumping 

case) (d 2 /dy 2 — a 2 )(d 2 /dy 2 — a 2 )4> = 0, a 2 = a 2 — iaRe. The boundary conditions are 

<py(±l) ± ra (±l)Kn = 0, 0(±1) = ±1. Similarly, with tp 2 = {D(y) + E(y)e i2a( - X ~^ + 

E*(y)e~ t2a ( x ~^} /2 , we have much more complicated equations (cf. Chu in [17]) and the 

boundary conditions 



D v (±l) + \[<t>y y {±l) + ^(±1)] = T Kn{i[0, ra (±l) + cf>* yyy (±l)} + D yy (±l)}, (4) 



E y {±l) + ro (±l)/2 = T Kn[0, ra (±l)/2 + £ ro (±l)], f?(±l) + 0,(±l)/4 = 0. After lengthy 
algebraic manipulations, we obtain = c e ay + cie~ ay + c 2 e ay + c 3 e~ ay , where c i: i = 0, 1,2,3 



are related to a, a, and Kn (cf. Chu in Ref. 17). 

To obtain a simple solution which relates to the mean transport so long as only terms of 
0(e 2 ) are concerned, we see that if every term in the x- momentum equation is averaged over 
an interval of time equal to the period of oscillation, we obtain 




where n is the integration constant and could be fixed indirectly in the coming equation 
below. Now, from Eqn. (4), we have 

£> w (±l) ± KnA^il) = -\[<t> yy {±l) + <j>; y (±l)] T Kn{i[^ ra (±l) + 0^(±1)]}, (6) 

where D y (y) = Hoy 2 + a±y + 02 + C(y), and together from the expression of ip2, we can 
obtain C(y) which also depends on the a, a, Re, Q,i = 0,1,2,3 (cf. Chu in Ref. 17). In 
realistic applications we must determine n from considerations of outlet conditions of the 
slab-region. a\ equals to zero because of the symmetry of boundary conditions. 
Once n is specified, our solution for the mean speed (u averaged over time) of super- 
flow is U = e 2 D y /2 = e 2 {C(y) - C{1) + R - KnC^l) + U [y 2 - (1 + 2Kn)]}/2 where R 
= — {[(f)yy(l) + 0*^(1)] — Kn[(j)yyy(l) + <f> * ( 1 ) ] } / 2 . To illustrate our results clearly, we adopt 
U(Y) = u(y) for the time-averaged results with y = Y in the following. 
Our numerical calculations confirm that the mean streamwise velocity distribution (aver- 
aged over time) due to the induced motion by the wavy elastic interface in the case of free 
(vacuum) pumping is dominated by i?o (or Kn) and the parabolic distribution — n (l — y 2 )- 
i?o which defines the boundary value of D y has its origin in the y-gradient of the first-order 
streamwise velocity distribution. Note that the Reynolds number here is based on the wave 
speed. The physical trend herein is also the same as those reported before 13 " 14 ' 17 for the 
slip-flow effects. The slip produces decoupling with the inertia of the wavy interface. 
Now, let us define a critical reflux condition as one for which the mean velocity U(Y) equals 
to zero at the center-line Y — 0. With the equation of U, we have n 0cr = Re(dp/dx)2 = 
-[a 2 Re 2 F(0)/200 + KnC'(l) - Rq]/{1 + 2Kn) which means the critical reflux condition is 
reached when n has above value. Pumping against a positive forcing greater than the crit- 
ical value would result in a backward transport (reflux) in the central region of the stream. 
This critical value depends on a, Re, and Kn. There will be no reflux if the pressure gradient 
is smaller than this Il . Thus, for some Il values less than n 0cr , the superflow will keep 
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moving forward. On the contrary, parts of the flow will move backward if n > n 0cr . 
As reported in Ref. 8, the rather small critical velocity (< 20/im/s) observed shows 
an apparent dissipation or attenuation of the superflow. Thus, we present some of the 
values of n (o;,-Re;Kn = 0,0.1) corresponding to freezed or zero-volume-flow-rate states 
(f\ U{Y)dY = 0) in Table 1 where the wave number (a) has the range between 0.10 and 
0.90; the Reynolds number (Re)— 0.1, 1, 10,50, 100. We observe that as Kn increases from 
zero to 0.1, the critical n decreases significantly (cf. Fig. 1). For the same Kn, once Re is 
larger than 10, critical reflux values n drop rapidly and the wave-modulation effect (due to 
a) appears. The latter observation might be interpreted as the strong coupling between the 
interface and the inertia of the streaming superflow. The illustration of the velocity fields 
for those zero-flux (zero-volume-flow-rate) or freezed states are shown in Figure 2. There 
are three wave numbers : a = 0.2, 0.5, 0.8. The Reynolds number is 50. Both no-slip and 
slip (Kn=0.1) cases are presented. The arrows for slip cases are schematic and represent the 
direction of positive and negative velocity fields. 

Some remarks could be made about these states : the transport being freezed in the time- 
averaged sense for specific dissipations (in terms of Reynolds number which is the ratio of 
wave-inertia and viscous shearing effects) and wave numbers (due to the wavy interface or 
vacuum fluctuations) for either no-slip and slip cases. Meanwhile, the time-averaged trans- 
port induced by the wavy interface is proportional to the square of the amplitude ratio 
(although the small amplitude waves being presumed), as can be seen in Eqn. (5), which is 
qualitatively the same as that presented in Ref. 11 for analogous interfacial problems. 
In brief summary, the entrained transport (either postive or negative and there is possibility 
: freezing) due to the wavy interface is mainly tuned by the n for fixed Re. Meanwhile, 
Ho cr depends strongly on the Knudsen number (Kn, a rarefaction measure) instead of Re 
or a. These results (cf. Table 1 and Fig. 1) might explain why there are rather small 
critical velocities for superflows in the temperature range where a supersolidity is observed 8 . 
Finally we like to stress that once the phase is solid then the elastically tensile stress should 
be larger than that of the liquid phase and the latter could be a crucial test for the existence 
of the supersolidity. 
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Fig. 1 Demonstration of Kn, Re and a effects on the IIo (zero-flux states). 
i?e is the Reynolds number (the ratio of the wave- inertia and viscous shearing dissipation). 
a is the wave number and Kn is the Knudsen number (a rarefaction measure). 



S 




Fig. 2 Demonstration of the zero-flux states : the mean velocity field U (Y) for 
wave numbers a = 0.2, 0.5, 0.8. The Reynolds number is 50. Kn is the rarefaction measure 
(the mean free path of the particles divided by the characteristic length) . 
The arrows are schematic and illustrate the directions of positive and negative U{Y). 
The integration of U(Y) w.r.t. Y for these velocity fields gives zero volume flow rate. 
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TABLE I: Zero-flux or freezed states values (Hq) for a flat wavy interface. 







Re 










Kn 


a 


0.1 


1 


10 


50 


100 





0.1 


4.5088 


4.5087 


4.5078 


4.4863 


4.4316 




0.2 


4.5269 


4.5269 


4.5231 


4.4496 


4.3275 




0.5 


4.6586 


4.6584 


4.6359 


4.4086 


4.2682 




0.8 


4.9238 


4.9234 


4.8708 


4.5714 


4.4488 




0.9 


5.0475 


5.0468 


4.9827 


4.6709 


4.5516 


0.1 


0.1 


2.3983 


2.3982 


2.3927 


2.2676 


1.9560 




0.2 


2.4003 


2.4000 


2.3774 


1.9532 


1.2217 




0.5 


2.4149 


2.4132 


2.2731 


0.7728 


-0.9054 




0.8 


2.4422 


2.4379 


2.0718 


-0.5885 


-3.4151 




0.9 


2.4542 


2.4486 


1.9825 


-1.1191 


-4.3972 
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